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Abstract 
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1 Introduction 



Boundary conditions have sophisticate influence on behavior of reaction-diffusion equa- 
tions. Since diffusion is the underlying mechanism for the spatial pattern formation in 
chemical reactions, spatial structure of solutions of reaction-diffusion equations can be 
sensitive to boundary conditions. Posing and realizing appropriate boundary condi- 
tions, e.g. for chemical reactions in open and large systems, is delicate (cf. Gray/Scott 
[7]). Furthermore, components of a system of reaction-diffusion equations can be im- 
posed with different boundary conditions. These make a big difference in behavior 
of systems and scalar equations, e.g. in Hopf bifurcations, spirals and other pattern 
formations. Typically multiple bifurcations occur more likely in systems. Moreover, 
stability of the bifurcating solution branches varies considerably from a scalar equation 
to a system. To distinguish influence of boundary conditions from that of interactions 
among the different species (components) of systems, we consider a scalar reaction- 
diffusion equation 

Au + Xu = f(u,X) in := (0,tt) x (0,tt). (1) 
We assume that the mapping / : R X R — > R is sufficiently smooth and 

/(0,A)=0, AJ(0,A) = 0, (2) 
i.e., it describes the nonlinearity of the problem (1) and implies that 

u = 0, A G R 

is a trivial solution of (1). 

We are interested in impact of symmetry-breaking in boundary conditions on bi- 
furcation scenarios. To this end we consider a square domain and impose the following 
conditions along its four sides 

Ou 

h o (fi)u(x,0) - h 1 (/i) — (x,0) = 0, 

% V 
ou 

h (/i)u(x, 7r) + hi(n) — (x, 7r) = 0, 

t , < 3 » 

=;«>,») = 0, 

Ou 

^,y) = o. 

Here h 0} hi : [0, 1] — > R are smooth functions and satisfy 

h o (0) = ^(1) = 0, h (fi) ^ for fie (0, 1], hii/i) ^0 for fi £ [0, 1). (4) 
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Figure 1: A homotopy between homogeneous Neumann and mixed boundary conditions. 

The boundary conditions (3) break the /^-symmetry into D 2 -symmetry. More 
precisely, properties of h (fi) } hi(fi) make (3) a homotopy from the homogeneous 
Neumann boundary conditions along the four sides of at ji = to the mixed bound- 
ary conditions at ji = 1, which are of the Neumann type along the sides x = 0, 7r; and 
of the Dirichlet type at y = 0, ir. 

Homotopy of boundary conditions has been used by Fiedler [3] and Gardner [4] 
to study global attractors and nonsingular solutions of a class of reaction-diffusion 
equations. They show that these solution sets are independent of boundary condi- 
tions. Nevertheless, if the equilibrium is nonhyperbolic and a bifurcation occurs, the 
bifurcation scenario, e.g. structure of attractors, may vary with respect to boundary 
conditions. This has been observed by Dillon/Maini/Othmer [2] in the study of pat- 
tern formation in generalized one-dimensional Turing systems and by Mei/Theil [13] 
in the analysis of steady state bifurcations as well as by Holder /Schaeffer [9] and Scha- 
effer/Golubitsky [14] on mode-jumping of von Karman equations. Using (3) as an 
example we study in this paper how reaction-diffusion equations react to a symmetry- 
breaking in boundary conditions. 

An outline of this paper is as follows. In Section 2 we consider variational form 
and symmetries of the equation (1). Section 3 describes bifurcation points of (1) along 
the trivial solution curve. In Section 4 the problem (1) at bifurcation points is reduced 
to algebraic equations via the well-known Liapunov-Schmidt method. We derive the 
bifurcation scenario at simple and double bifurcation points in Section 5 and illustrate 
these with a simple example in Section 6. 
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2 Variational Equations and Symmetries 

The classical regularity theory of elliptic problems ensures the C 2 ' s -H61der continuity 
of solutions of the linear problem 

Au + Xu = g 

with the boundary conditions (3) inside the square (cf. Wigley [15]). However, 
differentiability of solutions at the four corners depends strongly on properties of g. 
This linear problem is involved in the analysis of bifurcations of (1) with various right 
hand sides, here we want to avoid technicalities for the classical solutions and write the 
problem (1) into variational form. Namely, we study bifurcations of its weak solutions. 

2.1 Weak form 

We consider the Sobolev space 

X:=H\n) = {ueL>(ny, ^Yy eL2 ^} 

with the norm || • and for u,v £ X } ji £ [0,1) the parameter-dependent bilinear 
form 

bJu,v) := — f (ViiVi) + uv) dxdy — — [ [u(x, 0)v(x, 0) + u(x, n)v(x, it)] dx. (5) 

hi(fi) Jo 

For fx = 1 we choose the bilinear form 

bi(u,v) := — / (ViiVi) + uv) dxdy, (6) 

defined in the space X X X and 

X := {u £ i7 1 (0) ; u satisfies the boundary conditions (3) for fi = 1} 
The weak form of the linear problem 

Au — u = g in (7) 
with the boundary conditions (3) is 

Find ueH 1 ^), such that b^u, v) = (g, v) for all v £ H^tt). (8) 
In particular, weak form of the Neumann problem corresponds to ji = 0. 
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For domains £ C ' 1 , typically rectangle and L-domains (cf. Hackbusch [8, 
pp.118]), we have 

< C||u||i, n for all u £ H 1 ^), 



di 



dr 



< C\\u\\ 2 ,n for all u £ # 2 (fi), 



where C > is a constant. Note that 



< 



i = inf for all <j> £ H^dSl). 

H2(dn) U \ XI =4, 



Thus the bilinear form b„ 



is 



continuous and coercive on X X X. Moreover, if the 



inequality h (fj,)hi(fi) > holds, it is elliptic. Therefore, the problem (8) has a unique 
solution u(ji) £ X for every ji £ [0, 1) and all g £ H~ l (Jt). The solution u(ji) satisfies 
the boundary conditions (3). 

Denote the solution operator of (8) as 

TV) : g £ H~ l {ti) ^ T^)g = u(/i) £ H\Sl). (9) 

The operator TV) is linear and bounded for any fixed ji £ [0, 1). Furthermore, owing 
to the symmetry of Laplacian it is self-adjoint. This can be seen from the following 
equality for all /, g £ TT _1 (fi) 

(T(v)f, g)Hi{n)xH-i{n) 

{T{v)f)g dxdy 

V(T(fx)f)V(T(fi)g) + (T(fx)f)(T(fi)g) dxdy 

+ ¥r\r [(T(n)f)(T(n)9)\ o + (T( P )f)(T( P )g)\ ] dx 
ni[fi) Jo 1 y=o y=ir-i 

f (-A + I)(T(fx)f)(T(fi)g) dxdy + / (T((i)g)%-T((i)f ds 
h (fi) 



+ 



Jq [(T(rif)(T(rig) y=Q + (T(rif)(TMg) 
f ■ (T(f*)g) dxdy 



y=TT 



dx 



= if, T(/j)g) H -i(a)xm(a)- 
More importantly is the fact that together with the Gelfand relation 

the Riesz-Schauder theory is applicable to the operator TV) for all ji £ [0, 1), so that 
T(/i) has the following properties. 
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• Spectrum of the operator T(/i) consists of eigenvalues. There are maximally 
countable eigenvalues and can be ordered as 

Ai > A 2 > >0. 

The eigenspace associated to each eigenvalue A;, i = 1, 2, . . . is finite dimensional; 

• For i = 1,2,... the equation T(ji)u — Xu = f is solvable if and only if 
/ J_ Ker(TV) - A/). 

For fx = 1 we obtain the same conclusions with the bilinear form (6). In the sequel we 
consider the weak form 

G(u,\,n) ■=u + (X + l)T(n)u - T(fj,)f(u, A) = 0. (10) 

The mapping G : X X R — > X is obviously as smooth as / in (it, A). It is also 
continuously differentiable in ji due to the following properties of T(/i). 

Lemma 1 (Mei [12]) The operator T (ji) is continuous and differentiable with respect 
to fx in [0,1). Furthermore, the derivative u'(fi) = T\ji)g =: v(ji) for all g £ Y is 
given as the weak solution of 

Av — v = in 0, 
^(0,y) = 0, ^(7r,y) = 0. 

OX OX 

To calculate T(/i)'(/, we denote u = T(ji)g the weak solution of the equation 
Au — u = g with the boundary condition (3). Define 

Hw)'(-H" ( -'- (i2) 

It is easy to verify that v satisfies the boundary conditions in (11). Let v = w + v 
and substitute it into (11). We obtain the equation 

Aw — w = —(Av — i)) 

with the boundary conditions (3). Furthermore, 



< a "-"» = -Um 



-^ + y)g--u + 2(-^ + l) — 

IT IT TT Oy 
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Hence, the weak solution v of the equation (11) is 



-T(fi)(Av -v) + v 



Mi 1 ), 



2» 



-Tirfg + 2(^ - l)^-(T(n)g) + K - y] g 

IT IT oy \TT J 



(13) 



+ 



7T 



+ V T{n)g 



Remark: Note that T(/i) is self-adjoint, so is its derivative T\ji). Based on the 
formulation (11), one can calculate the higher order derivatives of u(ji) in a similar 
manner. 



2.2 Symmetries 

Let D4 be the dihedral group of the unit square and 

Si(x,y) = (1 - x,y), R(x,y) = (1 - y,x) 

be its generators. With Z 2 := {1,-1}, we define Z 2 X D4 = {iS; 8 £ _D 4 } and its 
actions on Y := L 2 (0) (d X) as 

ju(x,y) = ±u(^ _1 (x, y)) for all 7 = ±$, $ £ D4 and it £ Y. (14) 

The function spaces X, Y are obviously Z 2 X D 4 -invariant. Similarly, the L 2 -product 
is also Z 2 X D 4 -invariant. Corresponding to the boundary condition (3) we are partic- 
ularly interested in the subgroup 



D 2 :={S\,R 2 ; S\R 2 ,I}. 



(15) 



Let 



r := < 



' Z 2 X D4 if /(it, A) is an odd function in it and ji = 1; 
D 4 if /(it, A) is not odd in it and ji = 1; 

Z 2 X _D 2 if /(^ 5 A) is an odd function in it and /i ^ 1; 
, D 2 if /(it, A) is not odd in it and /i ^ 1. 

The T-equivariance of the mapping G, i.e., 

G(-fu } A, /i) = 7Cr(it, A, /i) for all 7 £ T, it £ X, A £ R. 

can be verified directly via the generators ±S\ and ±i?, respectively. 
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3 Bifurcation Points 



Since D x G(0,X,fji) = 1)^(0, A, /i) = and D u G(0, A,/i) = / + (A + l)T(fji) for all 
A, fx G R, a bifurcation occurs at a point (0,A,/i) on the trivial solution manifold 
{(0, A, /i); A G R, fi G [0, 1]} of (10) if the linearized problem 

D u G(0, \,n)u = u + (A + l)T(n)u = (16) 

has nontrivial solutions. By the definition (9) of T(/i) this equation is the weak form 
of the eigenvalue problem 

Au + \u = in = (0,tt) x (0,tt) 

with the boundary condition (3). To solve this problem with the rule of separating 
variables, we take the ansatz u(x } y) = Ui(x)u 2 (y) ^ and derive 

u" u" 

+ — + A = in = (0, tt) x (0, tt). 

Ui u 2 

Hence, u\ satisfies the equation 

u" + k\U\ = for some k\ G R (17) 
with the boundary conditions 14^(0) = u'^tt) = 0. Similarly, u 2 is a solution of 

u" + k 2 u 2 = for some k 2 G R (18) 

with the boundary conditions 

h o (fi)u 2 (0) — /ii(/i)u' 2 (0) = 0, 
h (fi)u 2 (ir) + /ii(/i)u' 2 (7r) = 0. 

These are eigenvalue problems of the one-dimensional differential operator with 
two different boundary conditions. Solutions of these problems are of the form 

(u\,ki) = (cos ra, n), n G N, 

(u 2 ,k 2 ) = ^ (/i) sin(A;(/i)y) + /ii(/i)A;(/i) cos(A;(/i)y), A;(/i)) , 
where A;(/i) G R satisfies 

2h (fi)h 1 (fi)k cos(kTi) + (hK/i) - hl{ji)k 2 ) sin(A;^) = 0. (19) 
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Thus eigenvalues of the Laplacian — A are given as 

\ = n 2 + k(fi) 2 (20) 

with the corresponding eigenfunction 

(f>(fi) := </>/||</>||, 4> '■= ( cos nx)[h (fi) s'm(k(fi)y) + hi(n)k(n) cos(k(fi)y)]. (21) 

As a function of the homotopy parameter /i, the function n(fJ>) has the following prop- 
erties (cf. Mei/Theil [13]). 

Lemma 2 Under the assumption (4), the equation (19) does not have integer solution 
for fi G (0, 1). Moreover, if [hi(fi) / h (fi)] < for all fi G (0, 1), then the solution k(fi) 
of (19) increases monotonously from m G N at ji = to (m | 1) £ N at /i = 1. 

Note that after multiplying the factor sin(A;7r) to the both sides of the equation 
(19) we can rewrite it as 

[h (fi) sin(A;7r) — hi(fi)k(l — cos(kir))][h (fj,) sin(A;7r) + hi(fi)k(l + cos(£;7r))] = 0. 

We use the following definition of parities of k(ji) in Ashwin/Mei [1], which are con- 
sistent with the parities of wavenumbers of the Neumann problem at ji = 0. 

Definition 3 The parities of the wavenumber k(ji) for Robin boundary conditions are 
defined as 

( EVEN if ho(fji) sm(k(fji)Tr) = h^k^l - cos(£;(/i)^)] ; 
k(h) = I for all fi G [0,1], (22) 

[ ODD if ho(j^) s'm(k(ji)ir) = —hi(ji)k(ji)[l + cos(k(ji)ir)] . 

In the rest of this paper we restrict the discussion to the case 

hi{jj) ~ ' \h {fi) J 
We conclude that bifurcation points of (10) on the trivial solution manifold are 

|(0, A(/i), /i); A(/i) = n 2 + k(/i) 2 , n G N, and k(/i) satisfying (19)}, (23) 

see Figure 2. The kernel Ker(_D„G(0, A(/i), ji)) is generically one-dimensional and 

Ker(_D„G(0, A(/i), /i)) = span[c/>(/i)]. 
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Figure 2: Bifurcation points of (10) in the parameter space (A, d). Here we have chosen 
h (/i) = /i, hi(jjL) = l-/i. 



From the statement (20) and the Figure 2 it is evident that two curves of bifur- 
cation points may intersect. In fact, these are generic as ji approaches zero and one, 
respectively. A intersection point corresponds to a multiple bifurcation of (10). At 
ji = 0, 1 solutions of the equation (19) are integers. In particular, at ji = the bound- 
ary conditions (3) reduce to homogeneous Neumann type. Thus the eigenvalues of — A 
are of the form 

A = n 2 + k 2 , n, k <G N U {0}. 
Such an eigenvalue is generically double and the associated eigenspace is 

Ker(_D„G'(0, A , 0)) = span[c/>i, c/> 2 ], 

where 

2 2 
(f>i = — cos(nx) cos(A;y), c/> 2 = — cos(A;a;) cos(ny). 

7T 7T 

Consequently, as ji — > 0, two different curves of simple bifurcation points approach the 
same double bifurcation point. However, with A as the bifurcation parameter, there is 
exactly one solution branch bifurcating from the trivial solution curve at every simple 
bifurcation point, while at a double bifurcation point for ji = the Neumann problem 
have up to four different nontrivial solution branches (cf. Mei [11]). On the other 
hand, as we have seen before, the _D 4 -symmetry of the Neumann problem breaks into 
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D 2 X _D 2 -symmetry as ji moves away from 0. This leads to the question how the 
bifurcation scenario reacts to such a symmetry-breaking in the boundary conditions. 
We want to investigate which solution branches of the Neumann problem and what 
symmetries of the solutions persist as the parameter ji varies. 

4 Liapunov-Schmidt Reduction 

To investigate solutions of (10) at a bifurcation scenario point (0, A , /i ) on the curve 
(0, A(/i), ji) in (23), we use the well-known Liapunov-Schmidt method to reduce 
the problem (10) to an algebraic system (cf. Golubitsky/Schaeffer [5] and Golubit- 
sky/Stewart/Schaeffer [6]). 

According to the Fredholm properties of D u G(0 } A , /i ) = / + (A + l)T(/i ), we 
have the decomposition 

X = Ker(D u G(0, A , /i )) © Im(D u G(0, A , /i ))- 

Suppose that the kernel Ker(_D„G(0, A , /i )) = span[c/>i, . . . , (f>i\ is ^-dimensional 
(£ = 1,2 generically). We write elements (u, A, ji) £ X X R X R as 

i 

u = + w = Z ■ (f> + w, 

8 = 1 

A = A + cr, 
jj = ii + u } 

where z = (z 1} . . . , z e ), <f> = (</>i, . . . , ^), z t} cr, v £ R and w £ lm(D u G(0 } A ,/i ))- 
Consider the projection Q := I — ^Cf =1 (c/>;, •)</>«• from X onto Im(D u G(0, A ,/i )). We 
rewrite the equation G(u } A,/i) = into a system 

QG(z ■ 4> + w, A + a, /i + v) = 0, (24) 
(I-Q)G(z-<f> + w, \ + <T, fio + v) = 0. (25) 

Solving w uniquely from (24) as a function of z, cr } v and substituting it into (25), we 
obtain the reduced bifurcation equation for z, cr } ji 

(I - Q)G(z ■ 4> + w(z, a, i/), A + a, fi + v) = 0. (26) 

We note that u;(0,0,0) = 0, D z w(0,0,0) = from (24). In the coordinate system 

z-(j>e Ker( J D„G(0, A , /i )) < — > z = (z u . . . , z e ) £ R^, 
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the operator equation (26) becomes a system of I algebraic equations 

B(z, cr, v) := l((f) t , G(z ■ cj> + w(z, cr, u), A + cr, /i + v)) . =0. (27) 

V / 2 = 1 

By definition the projection Q is T-equi variant. Thereafter the mapping B(z } cr } u) is 
also T-equivariant with respect to the induced action of T in R/, i.e., 

B(~fz, cr, v) = r yB(z, cr, v) for all (z, d,i/)eR'xRxR, 7 G T. 

We take the Taylor expansion of the components of (A + l)B(z, cr } v) = at the point 
(z, cr, v) = (0, 0, 0) and consider the truncated form 

-az t + (A + 1) 2 (^-, T'(fi )(z ■ $))v (28) 




= 0, i = l,...,L 

Here and in the sequel D uu f and D uuu f denote the derivatives of / at (it, A) = (0, A ); 
a G is a multi-index and 

Z l«l = 2 

represents the second order terms in the Taylor expansion of w at (z } cr } v) = (0, 0, 0). 

We recall that in the singularity theory a problem B = is finitely determined 
if there exists k G N, k < 00, such that the bifurcation scenario of B = and its 
k-jets jk(B) = 0, the Taylor expansion of B truncated at k-th order, are equivalent. 
The determinacy of a general problem at a bifurcation point is characterized by its 
reduced bifurcation equations. We refer to Golubitsky/Schaeffer [5] for more detailed 
discussions. For 3-determined problems solutions of the system (28) correspond one- 
to-one to those of the original problem (10) and contain all information of bifurcations 
of (10) at (0, A , fio)- 

With the knowledge of the bifurcation point (0, A ,/i ) and the kernel Ker(_D„Go), 
the terms (</>;, D uu f (z ■ c/>) 2 } and (</>;, D uuu f (z ■ c/>) 3 } in (28) can be calculated di- 
rectly. The other terms involve the derivatives T'(/i ) and D ZtZj w . Since the function 
w(z } cr } u) is defined implicitly by the equation (24), the term D ZtZj w is described as 
the unique solution v of the linear problem 

D u G v = QT(n )D uu f <t>i<t> h v G lm(D u G ). (29) 
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The term T'(fj, )(z ■ <f>) is calculated as a solution of the equation (11). In fact, via (13) 
we derive 

(<f>i, T'((i )(z ■<!>)) 

7T 



h(fi )((f> t 



+ y)T{fi ){z ■ 4>) 



+ T(fi )[lT(fi )(z ■ 4) + 2(^ - l)^(T(fi )(z ■ 4)) + (yl-y)( z . </>)]) 
h(fi ) 
h(fi ) 



-it ' 7T dy 
1 / , 2 , „ „ •_'// i) 



7T 



where h(fio) 



(Ao + 1) 2 



7T N 7T (7y 



5 Bifurcation Scenarios 

5.1 Simple bifurcations for /i G (0,1) 

Let A(/i) = n 2 + k(ji) 2 be a homotopy of simple eigenvalues of the Laplacian and 
A(0) = n 2 + ?7^ 2 , A(l) = n 2 + (ra + l) 2 - Then (0,A(/i),/i) is a curve of simple 
bifurcation points of (10) and 

Ker(_D„G)(0, A(/i), /i) = span[c/>] 

is 1-dimensional and (j> is given in (21). To obtain the generic bifurcation diagram of 
(10) at (0, A(/i ), Ho) for an arbitrary Ho £ (0,1), we consider the equation (28), i.e., 
3-jet of the reduced bifurcation equation, 

= -az + (\ + l) 2 (& T'([i )<t>)vz (30) 
+ (V, -D uu f (f> 2 ^ z 2 + (V, -D uu f (D zz w )(f> + -{D uuu f )(f> 3 ^ z 3 . 

Here z £ R is a scalar and 

2 



(</>, T'(fio)4) = (—) M 



(Ao + 1) 5 



7r ^ ^ 7T ^ dy^ 



Theorem 4 TTie problem (10) undergoes a pitchfork bifurcation at all points on the 
curve (0, A(/i),/i), ji £ (0, 1), i.e, the simple bifurcation points . Moreover, the truncated 
bifurcation equation (30) reduces to 



j 3 [(X + l)B(z,o-,u)] 



-a + av)z + cz" 



0, 



(31) 
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Wl 



th 



2h(ii ) 



TT 



2y 



TT 



d_ 

dy 



4> 



7;Duufo {D zz w ) (f> + -Duuufocj) 1 

2 



(32) 



Proof: It is easy to verify that the eigenfunction (j> of the Laplacian has the property 
((/>, (j) 2 } = 0. Thus the second order term in (30) vanishes. The conclusion follows 
directly from equations (30) and (31) consecutively. I 
The nontrivial solution of (31) is given as 

az/X 1 / 2 



G 



5.2 Double bifurcations of the Neumann problem 

For the Neumann problem (/i = 0) a generic double bifurcation point (0,A o ,0) has 
the property A = n 2 + A; 2 (0) with the wavenumbers n, k := A;(0) (n ^ k) as integers. 
Furthermore, we can choose 



Kei(D u G ) = span [(/>i, c/> 2 ] 



wi 



ith 



— cos(nx) cos(A;y), 



TT 



V2 



cos nx 



TT 



for n ■ k ^ 0; 
for n ^ 0, A; = 0; 
for n = 0, A; ^ 0, 



and 



</>2 := <f>i(y,x). 

On the other hand, taking into account the homotopy parameter ji in the boundary 
conditions, we see this double bifurcation point is split into two simple bifurcation 
points (0, A 8 -(/i), /i), i = 1,2 with Ai = n 2 + k 2 (ji) and A 2 = A; 2 + ^ 2 (/") for /i / 0. We 
want to investigate bifurcation scenario of (10) at a double bifurcation point (0, A , 0) 
and its variation with respect to the homotopy parameter ji. 

Note that (</>;, </>j</>;) = for all i,j,l = 1,2. The equation (29) can be solved 
analytically (cf. Mei [11]). Together with the statements 
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Att 2 

3 
2^ 



for n ■ k ^ 0, 
for n ■ k = 0, n 2 + k 2 ^ 
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and ((1)1 $) 



IT' 



for i ^ j, we simplify the equations (28) into 



-az 1 + (A + 1) 2 (</>i, T' (no)(zi(j)i + z 2 $ 2 ))v + c x z\ + c 2 z x z\ = 0, 
-az 2 + (A + l) 2 (</>2, T'(/i )(^i(/>i + z 2 cf) 2 ))u + c 2 z 2 z 2 + c x z\ = 0. 

Here ci, c 2 are constants. More precisely, if n ■ k ^ 0, we have 

1 



Cl 



c 2 



6ir 2 
1 



9 1, , 2 45(P-n 2 ) 2 + 4Pn 2 

4 uuu/ ° 4 l " u/0j (F - 3n 2 )(n 2 - 3F)(n 2 + F) 



6ir 2 

If n = 0, k ^ 0, then 



3D /o - Q(DuufoY 



1 



(P -n 2 ) 2 -4Pn 2 



n 2 + k 2 \[(k 2 + n 2 ) 2 — 16k 2 n 2 ] 



V 



1 



Cl 

c 2 



6ir 2 
1 

2^ 



3 5 2 

-D uuu fo + -^-^(Duufo) 



(33) 



If i/ = 0, the equations (33) coincide with those in Mei [11], and yield four nontrivial 
solutions of (10) with symmetries of the isotropy subgroups of </>i, (f> 2 and (f)\ ± c/> 2 , re- 
spectively. For v ^ 0, the forced symmetry-breaking in boundary conditions introduce 
in (33) the terms 

(A + l) 2 (^, f(/i )Mi + ^))i' 



(0) 



2 (-y 



7T 

1,2. 



<9y 



Here, 



rfi = </ 2 = -(— )(0), ifre-fc^O, 

cfi = — C— ")(0), d 2 = 0, ifn = 0, k ^ 0. 
it \g x ' 



The system (33) reduces to 

[ - a + d]_v + c x zf + c 2 zf ] ^i = 0, 

[ — a + o? 2 i/ + c 2 z\ + c 1 z^z 2 = 0. 



(34) 
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Remark: The coefficients in the equations (32) and (34) are related as follows have 
properties 



lim a = di for </> = </>;, i = 1 , 2 

/i-»0 

lime = c\. 

At->0 

Solutions of the system (34) are 

- diV\l/2 \ / /(7 — d 2 V\l/2 




(ci - c 2 )a + (c 2 o?i - c 1 d 2 )i / \i/2 



(35) 



These lead to four bifurcating solution branches of the original problem (10), i.e., (1). 
The solutions in (35a) are pure mode solution branches with the isotropy groups of 
c/>i, c/>2, respectively. They correspond to those bifurcating at the simple bifurcation 
points on the curves (0, A 8 (/i), /i), z = 1,2. The solutions in (35b) involve both (f)\ and 
(f>2 modes. They are called the mixed mode branches. 

If n ■ k = 0, the terms d\v, d 2 v break the D 4 -symmetry of the Neumann problem 
and the mixed mode solution branches have merely the trivial symmetry. Moreover, 
the pure mode and mixed mode solution branches may intersect at 

(ciO?i — c 2 d 2 )v (c 2 o?i — c x d 2 )v 

a = or a = , 

(ci-c 2 ) (ci - c 2 ) 

and induce a secondary bifurcation, respectively. 

If n ■ k ^ 0, then d\ = d 2 and the _D 4 -symmetry is preserved in (34) and the mixed 
mode solutions (35b) becomes 

^ ( fU — 0?iI/\l/2 ^ jU — diV\l/2 



Ci + C 2 ' V Ci + c 2 

Symmetries of these solutions are the isotropy groups of the eigenfunctions (f)\ + (f> 2 
and (f)\ — 4> 2} respectively. 

As an conclusion, we see that all four bifurcating solution branches of the Neumann 
problem at a double bifurcation point persist if we vary both A and ji as bifurcation 
parameters. Moreover, symmetry of these bifurcating solution branches is preserved 
for those with the wavenumbers n ■ k ^ and is broken for those with the wavenumbers 
n = or k = 0. 
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6 A Simple Example 



Choose 

f(u,X) = X(u 2 + u 3 ). (36) 

We consider bifurcation scenarios at the corank-2 bifurcation points u = 0, /io = 
and A = 5, 10, 20, respectively. Moreover, we examine variations of the bifurcation 
scenarios as homotopy parameter ji moves away from zero, i.e., as the homogeneous 
Neumann boundary conditions with D 4 -symmetry are perturbed. To simplify the 
discussion, we take h (fi) = /i, hi(fi) = 1 — ji. Note that at ji = we have 

Kei(D u G ) = span [(/>i, c/> 2 ]. 

and the inequalities C\ ^ 0, c 2 ^ and c\ — c\ ^ holds for all n, k G N U {0}. 
(1) Wavenumbers n = 1, k = 2: For A(0) = 5 and n = 1, k = 2 we have 

2 2 
</>i = — cos x cos 2y, (/) 2 = — cos 2x cos y. 

7T 7T 
T, 1 , 4 5695 110220 o i • / X • /orx i 

furthermore, a = —, c\ = -, c 2 = — . Solutions (zi, z 2 ) m (35) becomes 

7T 1327T 2 1327T 2 

/ /1327r(7rcr - 4i/x 1/2 \ / /1327r(7rcr - 4i/)^i/2^ 

l ± ' 5695 J ' °J' 1°' ± ^ 5695 



/ j 



1327T(7TCr - 4l/)x 1/2 /1327T(7TCr - 4l/)x 1/2 



± U 115915 ' ' ± V 115915 



Figures 3 and 4 show the pure and mixed modes solution branches. 

(2) Wavenumbers n = 0, k = 1: We have A(0) = 1 for n = 0, k = 1, moreover, 

^ ^ x ^ 

(pi = cos y, q>2 = cos x. 

IT IT 

4 19 3 

Simple calculations show d\ = — , d 2 = c\ = — -, and c 2 = — . The solution 

7T 6lT 2 IT 2 

branches described by (zi, z 2 ) in (35) are 

/ /Qtt (ira — Au)\ 1/2 \ / /67r 2 <7\i/2 N 

±(— ^ L ) ' , 0, ± ' W 



19 ' ) V V 19 

'67r(7T(T — 76l/)\ 1/2 /67r(7T(T + 72l/)\l/2 N 



37 ' 1 V 37 

The pure </> 2 -mode solution branch meets a mixed mode solution at a = 76v/ir and 
induces a secondary bifurcation. 
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Two pure mode solution branches 




Figure 3: Two solution branches of pure (f)\ and (f> 2 modes. 




Figure 4: Two solution branches of mixed (f)\ and (f> 2 modes. 
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